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Abstract. - Phononless plasmon assisted transport through a long disordered array of fi- 
nite length quantum wires is investigated analytically. Generically strong electron plasmon 
interaction in quantum wires results in a qualitative change of the temperature dependence 
of thermally activated resistance in comparison to phonon assisted transport. At high tem- 
peratures, the thermally activated resistance is determined by the Luttinger liquid interaction 
parameter of the wires. 



Hopping transport in a quasi-one-dimensional system formed by a parallel arrangement of 
conducting wires is of much relevance to a number of experimental setups, including quantum 
wire arrays in hetero junctions [1], carbon nanotubc films [2], atomic wires on silicon surface [3], 
and stripe phases [4]. At finite length of constituent wires, such systems represent particular 
examples of granular arrays, where a one-dimensional wire plays the role of a grain. Considered 
as a granular array, the array of parallel quantum wires is rather peculiar because of the very 
long charge relaxation time in a one-dimensional wire. Due to this peculiarity, the theoretical 
description of thermally activated transport in arrays of long quantum wires requires taking 
into account the charge dynamics and treatment of the interactions beyond the capacitive 
model adopted in recent theoretical investigations of transport through disordered granular 
arrays [5]. 

In this letter we show that charge-density fluctuations (plasmons) in the array can act as 
the agent promoting thermally activated transport, thus providing the possibility for phonon- 
less inelastic transport. As the result of generically strong plasmon-electron coupling in a 
quantum wire, the features of plasmon and phonon assisted transport are qualitatively differ- 
ent. We provide a qualitative explanation of plasmon assisted transport, identify the transport 
regimes, where the features of plasmon and phonon assisted transport are either similar or 
substantially different, and derive analytic expressions for the temperature dependence of the 
thermally activated resistance for a special model of a strongly correlated disordered array of 
quantum wires. 

The model we formulate below is special, because it combines two seemingly incompati- 
ble features: i) it is strongly disordered for single electron transport; ii) it is much weaker 
disordered for propagation of plasmons. 
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Consider a one dimensional array of parallel identical quantum wires of length L and 
diameter a placed regularly with the interwire distance d, L ^ d a (see Fig. We 
investigate transport in the direction perpendicular to the wires. The spectrum of low-energy 
plasmons in a single isolated wire is equidistant with energies Ei,n = '^n, where L is the 
length of the wire, and Vi is the plasmon velocity along the wire i. For identical wires at regular 
positions, the intra- and interwire interactions between the charge density fluctuations do not 
change along the array. Then each localized plasmon level broadens into a plasmon band with 
truly continuous spectrum, quite analogously to the formation of electronic bands in the tight 
binding model. The role of the hopping in the tight-binding model for plasmons is played 
by the matrix element of the charge-density interactions between the neighbor wires. The 
formation of plasmon bands is reflected by the dependence of the plasmon velocity along each 
wire on the wave vector p along the array (that is, perpendicular to the wires) vq u{p). The 
particular form of plasmon dispersion depends on details of the interwire interactions, yet the 
function u{p) should be periodic with a period of one Brillouin zone. That is why we choose 
the specific form 

Up ~ vq + vi cos(7rp), — 1 < p < 1. (1) 

The chosen form of dispersion can be considered as the first two terms of the Fourier expansion 
of some general dispersion law. The plasmon energy within a band centered around the level 
n is given by e„(p) = 

Furthermore, if the potential energy barriers between the wires are high, the electronic 
bands for the motion along the array (in Y-direction in Fig. are very narrow with the 
width t w 2m {a\+i) where A describes the decay of the single particle wave function 

inside the barrier, is the effective mass of an electron. Now let us introduce the disorder as 
a random height of the energy barriers between the neighbor wires. Such a disorder induces 
fluctuations of the decay parameter A thus rendering t random. As the result, the single 
particle wave functions become localized. At (SX^) « ^g^^"^^'', the localization length of the 
single particle wave function is of order d, that is, an electron is localized within a single wire. 

The resistance of a long disordered one-dimensional array is determined by so-called breaks, 
the junctions between two neighboring wires with exponentially high resistance [7]. Let us 
denote the energy cost to transfer an electron over the break as Ea- To facilitate the trans- 
port over the break, the energy Ea should be borrowed by absorption of a bosonic excitation. 
Now we show that due to the transfer of energy between charge density excitations in dif- 
ferent wires, a bosonic bath with continuous energy spectrum can be created at the break. 
Then any energy deficit Ea can be matched by absorption of a plasmonic excitation, which 
constitutes the microscopic mechanism of plasmon assisted transport [8]. For two isolated 
wires forming the break, the matching condition Ei_n — Ea cannot be satisfied for arbitrary 
Ea because of discreteness of Ei,n. However, due to the charge-density interwire interac- 
tion Hint = Si X^n n' whcrc n,n' mark the plasmon modes, the energy 
can be transferred between excitations localized in different wires. Treating the interwire 
interaction perturbatively, we deduce, that the transfer of plasmon energy from the nearest 
neighbor of the wire i leads to splitting of the plasmon level _Ej „ into M sublevels with ener- 

gies e- „ = Ei^n — e-^]^]^e- ' where I = 1,M runs over the plasmon levels of the wire i + 1, 
and M denotes the number of levels that have a nonvanishing interaction strength with the 
level {i,n). Each split energy sublevel corresponds to a complex plasmon excitation, which is 
a superposition of plasmons localized in the two wires. Note that all the split sublevels are 

|y|2 

situated in the energy window ^-j^ around i?i^„, where 1^ is a typical value of the interwire 
interaction, and 6i is the plasmon inter level distance in a single wire, \V\ <C i5i. Now we take 
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into account that each energy level of the wire i + 1 is in turn split due to the interaction with 
the next-nearest wire i + 2. Repeating the previous calculation with -Bi+i,; substituted by the 
energy of a split level ^ we conclude that the energy interval |yp/i5i around Ei^n is now 
filled with levels. By recursion we obtain that the transfer of energy from n neighbors 
of the break results in M" plasmon sublevels within the energy interval around 
Therefore, due to the transfer of energy from plasmons localized in distant wires to the break, 
a bath of bosonic excitations at the break is formed that allows to match the energy Ea with 
the higher accuracy the more localized plasmons are involved in the creation of a bosonic 
excitation. 

In the disordered array the number of wires that can supply energy to the break is limited 
by the plasmon localization length that can be evaluated as = UgTp, where Tp denotes 
the plasmon mean free time and Ug is the plasmon group velocity. The mean free time Tp is 
related to the fluctuations of the matrix element of the interwire charge density interactions 
{{SVi^i+i)"^) — (2^)d^°t/ T d ' where z/p = ^-^^^^ is the plasmon density of states in the middle 
of the band n, and lu is the plasmon frequency. Furthermore, the random height of interwire 
tunneling barriers affects the interwire charge density interactions ~ e^a^|-!/'iP|V'i+iP/rf 

through the randomness of the wave functions ipi. Another source of randomness in the 
model is a random charging energy due to the random environment (local concentration of 
charged impurities) around each wire. Random energy barriers are directly felt only by the 
exponentially decaying tail of the single particle wave function within the barrier, whereas 
random fluctuations of the charge density interactions are mostly determined by the part of 
the single particle wave function within the wire. That is the reason why the disorder that leads 
to the localization of single particle wave functions is felt only weakly by the charge density 
fluctuations. Using the normalization condition (a + e~^^'^dx) = 1 one can relate 
the fluctuations of the decay factor Xi with the fluctuations of | -0^1^ that in turn determines the 
fluctuations of the charge density interactions and the plasmon mean free time. One obtains 
Tp = leT^^^i^fl^a^^^A^) • regime, when electrons are localized on the length d, we obtain 

Ip/d = wg(aA + l)''e^^'^/ {2LuXd)^ . For Ad > log ( .„ ^al+i)'^ \ plasmon locahzation length 
greatly exceeds the single particle localization length. If the length of array is less than Ip, 
the plasmons propagate freely along the array, whereas the electrons are strongly localized. 
In that regime, the dispersion can be used to describe the propagation of plasmons. 

Treating the interwire charge density interaction perturbatively, we can write the transition 
rate with the absorption of a plasmon similarly to a transition with the absorption of a phonon 
using the Fermi golden rule [9] 

/CO 
'^pY. Y1 \Vn{p)\NB{en{p))n-E,n)[l- f{Ea + Ek)]5{e,,[p)-Ea-E„,-Ek). 

n m,k=0 

(2) 

Here Vn{p) is the strength of interwire charge density interaction for the plasmon mode n, 
NBi^nip)) is the occupation number of the plasmon mode, f{E„i) denotes the Fermi dis- 
tribution and describes the occupation of the m-th single-particle energy level in the wire, 
Em = ^^m. For narrow plasmon energy bands, the perturbative approach suggests that if 
the energy Ea lies in the gap between the plasmon bands the hopping over the break is blocked. 
This suggestion turns out to be wrong because of a conceptual difference between the plas- 
mon and phonon transport mechanisms. Whereas the phonons represent a bath of bosonic 
excitations that is independent of electrons, the plasmons are "made" of electrons themselves. 
Consequently, while the electron-phonon interaction can generally be treated perturbatively, 
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the perturbative treatment of plasmons is possible only under special conditions. The appli- 
cability of the perturbative treatment of plasmons is determined by the relation of two time 
scales: the characteristic time of plasmon dynamics tp and the characteristic time of a single 
electronic hop th ■ If the Coulomb interaction in a grain is well-screened or the plasmons are 
strongly localized, then tp is the characteristic relaxation time of a plasmon excitation within 
a single grain. For tp ^ th the plasmons can be neglected in transport. The description of 
interactions thus reduces to the capacitive model [5] . For the delocalized undamped plasmons, 
the time tp is associated with the formation of an extended in space plasmonic excitation. In 
that case, tp ^ th correspond to the regime of a strongly nonlinear coupling between plas- 
mons and electrons, and the perturbative treatment of plasmons is incorrect. Plasmons in 
one-dimensional wires represent a profound example for that regime. For the model c;onsid- 
ered in this paper tp ~ L/vi, and at the break the condition tp -C th is always satisfied. 
As we show below, due to the strong electron-plasmon coupling, the nonlinear effects lead to 
the creation of plasmon complexes with energies covering the whole spectrum continuously, 
even though the plasmon bands initially are very narrow. This in turn leads to plasmon as- 
sisted transport with a temperature dependence that is qualitatively different from the case 
of phonon assisted transport. In the regime tp ^ th, the effective interaction time is limited 
by th- Then the plasmon dynamics is essentially independent of the electron dynamics and, in 
the case of a continuous plasmon spectrum, the plasmon assisted transport is quite analogous 
to the phonon assisted one. 

A unique feature of the chosen model is the applicability of the bosonized description 
that allows exact treatment of interactions and hence nonperturbative treatment of plasmons. 
Precisely, the plasmon dynamics is described by the action 

S= dp dT —I -IdrBpl^ + up\d,ep\^ \ , (3) 



1 Jo J-L/2 t Up 

representing a finite size generalization of the sliding Luttingcr liquid model [10]. The relation 
of the plasmon velocity Up and the Luttinger liquid constant Kp with inter- and intrawire 
interactions has been calculated in [10]. A fermion annihilation operator in the wire n, ^n{x), 
is represented as 



dp(t)^{x)e 



— tTT'pn 



(4) 



where % = i?, i denotes the chirality, (f>^{x) is a chiral bosonic field, and is a Klein 
factor. The chiral field (j)^ is in turn expressed through the field Qp{x) and its dual $p(x), 

The resistance of the array is calculated along the lines of Rcf. [7]. Let us parameterize 
the tunneling matrix element between two wires in the form = exp(— The 

parameter y can be associated with an effective distance between the two wires. This effective 
distance is random, its distribution follows from the distribution of the heights of potential 
barriers. Since a break, being a junction with exponentially large resistance, is not shorted 
by other resistances connected in parallel, we can write the resistance of a break in the form 
Ri = Rq cxp[2|?yi.i_|-i T)]. Here Ea denotes an addition energy to transfer an electron 

over the break. We remind that the disorder enters the model only as a random distribution 
of addition energies Ea and effective distances The function f{Ea,T) accounts for the 

effect of thermally activated plasmons in the resistance of the break. According to Ref. [7] , 
the probability density p{u) for the resistance R/Rq = e" is proportional to e~^"^, where A is 
the area in the {y,Ea) phase space that results in the resistance e", and g is the linear density 
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of localized single particle states. The resistance is calculated a.s R — R^ly du e"~S"^("\ 
where ly is the length of the array. Therefore, in order to calculate the resistance of the array 
in the localized regime, we have to obtain an expression for the resistance of the break 
Since the break is not shorted by other resistances, we conclude Ri = l/ci, where cti is the 
conductance of a break. Assume that the break is formed by a junction between the wires 
with numbers and 1. We take the position of a pinhole connecting the two wires as x. 
In the linear response approximation the current through the break Ii is determined by the 
correlation function [11] 

X{t) = (Tr (vl/o(a;,r)vI/t(a:,r)vI.i(x,0)*S(x,0))) (5) 

that characterizes the probability of a single hop over the break, X+{t) — X{t > Q) and 
X-{t) = X{t < 0). Here toi is the tunneling matrix element. In the bosonized representation 
(0J, the correlation function X{t) factorizes in the correlator of Klein factors and the correla- 
tor of bosonic exponents that we denote as Xf,(r). The time dependence of the Klein factors 
F„_^(r) is given by the ground state energy of the wire n that includes the capacitive inter- 
action between the wires. Thus the correlator of the Klein factors is proportional to e"^""^. 
Furthermore, calculating the correlators of bosonic exponents, we cast the expression for the 
current into the form I{eV) oc du)' Jo{eV +uj')Z{—uj'). Here Jo(w) ~ 5{uj — Ea) — 5{uj + Ea) 
defines the zero-temperature current, whereas the influence of thermally activated plasmons 
is contained in the factor Z{ijj) that is given by the the Fourier transform of 



Z{t) = exp 



CT— ±1 m— 1 



(6) 



Here Kp ~ Kp + 1/Kp, and the average over the plasmon wave vector p is defined by 
= jj^ • (1 — 2 cos(7rp) -I- cos(27rp)) dp. Further we assume the coupling constant to be 
p-independent, Kp = n and use the simplified dispersion law JQ). Note, that in approximate 
evaluations it is much more important to keep the p-dependence of the velocity Up that re- 
flects the formation of plasmon bands than the p-dependence of the coupling constant Kp. 
To the lowest order in p-dependent terms, the latter just leads to the averaging of the single 
Luttinger liquid result over the coupling constant. The basic average to be used in subsequent 
calculations reads 

loibvi) + {I - j^)Ii{bvi) , 



bvi ' 



(7) 



where b = (m/3 ± 11)11/ L, and Ii,{z) denotes the Bessel function of complex argument. For 
large times t, gives asymptotically 



e 



where w ~ vq — v\. Despite being obtained for v\ < wq, 1® is essentially nonperturbative 
in v\. The relevant values of the transport time are restricted by the hopping time t^. For 
v\\mfi ^ iThl < 1 the correlations giving rise to plasmon bands do not develop, and the short 
time expansion of has to be used instead of ^ . The latter is equivalent to the perturbative 
treatment of plasmons, leading to a result similar to the phonon mechanism of hopping. 

At low temperatures, T <C the major contribution to Z{t) in (|SJ) is given by the term 
with m — 1. Leaving only that term, expanding the logarithms, and substituting ^ in © 
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we obtain 



oo 




) 



■+1 




e-/3(^(«+0+l' 



sm 



(9) 



2 



where ei„(aj) — ^{l — n), and ly 



I, n for LO — ein ^ respectively. Each term in © describes a 



thermal excitation of a multiparticle plasmon complex with a continuous density of states. The 
plasmon complexes described mathematically by © form the bath of bosonic excitations that 
facilitate transport over the break. Higher values of to in @ would correspond to excitations 
involving progressively more plasmon modes. The leading terms in at low temperatures 
are given hy l^n — 0, 1. Those terms result in the leading low-temperature contribution to tri 
in the form 



Further calculations closely follow Ref. [7] leading to the temperature dependence of the 
resistance 



Since in the low-temperature regime the temperature broadening of discrete plasmon levels 
in a single wire is less than the interlevel separation, the finiteness of the wires essentially 
determines the physics of transport. In particular, the LL interaction parameter k does not 
enter the temperature dependence of the resistance, except as a general prefactor in 
For comparison, in the case of phonon-assisted transport in the same model, as well as by a 
perturbative treatment of the plasmon-electron interactions the preexponential factor in 
goes like T^/^. The result Hll|l also differs from the thermally activated resistance of a single 
disordered wire obtained recently in Ref. [6] . 

At high temperatures, when T 3> the temperature broadening exceeds the interlevel 
separation in a single wire. Then the plasmons in the array behave similarly to plasmons in an 
array of infinite length wires. The conductance of a single junction depends on the addition 
energy not exponentially but as a power law ai '--j Ea ■ [max (T, Therefore, the 

resistance of the array is not determined by the resistance of a single break, it is rather given 
by the average over the resistances of all junctions. The fact that the leading contribution to 
the conductance at high temperatures is given by the result for infinitely long wires implies 
that the coherence of the single particle motion is broken already by a single hop between 
the two neighbor wires. The latter justifies the averaging over the all junctions in calculation 
of the resistance. The expression for the resistance can be written as a Drude formula with 
the interaction and temperature dependent mean free time Tf. The expression for the mean 
free time can be organized as an expansion in powers of the small parameter e . At a 
typical charging energy Ea > T, the leading term in the expression for the mean free time 
is temperature independent, r/ oc 1/{E'^~'^). For comparison, phonon assisted transport in 
that temperature regime still has a thermally activated character with the preexponential 
factor T~^/-^ in At temperatures even larger than the typical charging energy, Ea < 

T, Tf exhibits the power-low temperature dependence typical for transport across a sliding 
Luttinger liquid, r/ oc T^~'^ / (Ea) [10]. Therefore, at high temperatures, the mean free time 
is determined by the Luttinger liquid interaction parameter k. The calculation of the phonon 
assisted resistance in that regime using the Fermi golden rule leads to the Drude formula with 
the logarithmic temperature dependence of the mean free time Tf oc InT. 




(10) 




(11) 
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Fig. 1 - Geometry of the model. The arrow shows the direction of the current. 

In conclusion, we demonstrated the possibility of plasmon assisted inelastic transport in 
the particular case of a disordered granular array, an array of parallel quantum wires. Due 
to large charge relaxation time in a wire, the plasmon-electron interaction has to be treated 
nonperturbatively. In the result, thermally activated resistance has qualitatively different tem- 
perature dependence for plasmon assisted transport as compared to phonon assisted transport. 
Despite the specific quasi one-dimensional geometry of the grains in this model, the described 
phenomenon is believed to be of general importance for granular arrays with delocalized or 
weakly localized plasmons. 

* * * 
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